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Abstract.- We study the influence of the different choice of unit cells on the Bloch solutions 
of Schrodinger equation for one-dimensional periodic Kronig-Penney models with rectangular 
potential barriers or potential wells and partially constant effective mass. We obtain generalized 
Kronig-Penney relations for bulk states and exact expressions for the corresponding periodic 
parts of the Block wave functions for the two possible choices of the matching conditions in 
an unit cell. We show also that our analytic expressions reduce to the well known expressions 
for the Kronig-Penney relations and Bloch waves under constant effective mass and appropriate 
matching conditions 
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The superlattice (SL) electronic structure is of primary importance and interest in both 
fundamental and applied aspects, and it has been extensively studied recently [1-10]. Theoretical 
studies of the electronic properties of SLs may be carried out in the framework of one-dimensional 
Kronig- Penney (KP) models [4,11], in particular, calculating the electronic surface states [2,4- 
6], the electronic interface states [3,7-9] and the electronic bulk states [10] in the semiconductor 
superlattices. The basic point in these considerations is to determine the general solution to the 
corresponding one-dimensional Schrodinger equation together with the conditions for existence 
of surface and interface states, which conditions depend strongly on the corresponding wave 
functions [2-9]. It hhas been noted by several authors [12-14] that the kind of the wave function 
depends on the choice of the unit cell, therefore, it seems very important to have a more detailed 
knowledge of this dependence between the choice of the unit cell and the wave function. 

The aim of this letter is to consider some aspects of this problem, in particular, how the 
periodic parts of the Bloch solutions to the corresponding Schrodinger equation with periodic 
potential and partially constant effective mass depend on the choice of the unit cell, more 
precizely, on the matching points for the wave functions and their derivatives in this unit cell. 

We are going to consider two cases of infinite one-dimensional KP-models: with rectangular 
potential barriers given on Fig.l, and with rectangular potential wells given on Fig.2. In the 
first case (potential barriers) the potential V{x), determined by the conditions: 

iv = const, n^{x) : nb+ {n + l)a < X < {n + l){a + b), 

[O, n'^{x):n{a + b)<x<nb+{n + l)a, (n = 0, ±1, ±2, . . . ) 

is periodic with respect to x with a period equal to the lattice constant (a + 6), i.e. V{x) = 
V{x + (a + 6)). For an electron with a partially constant efi^ective mass m{x) 

I mi = const, Q.^(x) : nb + (n + l)a < x < (n + l)(a + b) 
m[x) = S 9 (2) 

I 7712 = const, fl (x) : n{a + b) < x < nb+ {n + l)a, (n = 0, ±1, ±2, . . . ) 

moving in the periodic potential field (1), the one-dimensional stationary Schrodinger equation 

(f^(x) 2m(x) / X 

has, according to the Bloch-Floquet theorem [4,11], the following periodic in x solutions with 
period {a + b): 

ipi{x) = A exp [{—ia + j)x] + B exp [{—ia — 7)x] (4) 

in "regions I" : r2^(x) : nb + {n + l)a < x < {n + l)(a -|- b), and of the form 

ipiiix) = C exp [i{—a + P)x] + D exp [i{—a — P)x] (5) 

in "regions 11": Q,'^{x) : n{a + b) < x < nb + {n + l)a. Further in our consideration we assume 
the quantities (3 and 7, introduced by (4) and (5), to be real and positive and their explicit 
dependences on the electron energy E,{{) < E <V), are given by 



(6) 
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and a is determined by the matching conditions. To obtain the coefficients A, B, C and D of the 
wave functions (4) and (5) we make use of the matching conditions at three successive points 
(xi,X2,X3) in one unit cell such, that \xi — x-^\ = (a + b), and \xi — = a, \x2 — x^l = b, or 
\xi — X2\ = b, \x2 — x^l = a, and we assume xi < X2 < X3 (Fig.l). Since there exist only two 
ways to choose these three points on the x axis we introduce them by the following relations: 

xi + a = X2, X2 + b = xs, (7) 
xi + b = X2, X2 + a = xs, (8) 

Relations (7) and (8) specify the two types of unit cells: the unit cell of type "well-barrier" 
which we call KP-1 type cell (Fig. la), and the unit cell of type " barrier- well" , called KP-2 
type cell (Fig.lb), respectively. Thus, if we use the matching conditions used in [1,3,15-17] and 
denote 2/ = ^ > the solutions (4) and (5) for the unit cells (7) and (8) must obey the following 
conditions aX, x = X2- 

1pl{x2) = 1pIliX2) 

dipi dipii (9) 

= 

It follows from the periodicity of il^{x) that the solutions (4) and (5) and their first derivatives 
must satisfy also at the points x = xi and x = x^ the conditions 

di^i. . di,ii, . (10) 

for the KP-1 type unit cell defined by (7), and the boundary conditions 

'^i{x\) = ipnixs) 
dipi dijjii (11) 

for the KP-2 type unit cell defined by (8). Prom (9) and (10) and from (9) and (11) we obtain 
two systems of four linear homogeneous equations each, for the four constants A,B,C and D, 
reflecting the differences of the two types unit cells. These two systems have nontrivial solutions 
only if the corresponding determinants are equal to zero, which leads to the same equation for 
/? and 7 given by 

cos [a{a + b)] = ch{jb)cos{pa) + ^ ^ ^ —sh{jb)sin{(3a). (12) 

This equation (12) we call generalized Kronig-Penney relation. 

Since only three of the four coefficients A, B, C and D are linearly independent we can 
express any three of them through the fourth. So, when A = fi-'^B and C = ^''jjD, where the 
index j = 1,2 denotes the two types unit cells, the corresponding periodic wave functions (4) 
and (5) become 

il)^j{x) = B ^fx^jexp[{—ia + j)x] + exp[{-ia — j)x]^ , (j = l,2), (13) 
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ilPjjix) = D exp [i{—a + I3)x\ + exp a — , (j = 1, 2). 



(14) 



For the KP-1 unit cell the corresponding nonhomogeneous system of equations (9) and (10) 
gives the following expression for 



cos{Pa) — 




sin{Pa) — exp [—ia{a + b) — 76] 


—cos{(3a) + 




sin{(3a) + exp [—ia{a + 6) + 76] 



exp(— 27x2) 



(15) 



Similarly, the nonhomogeneous system of equations (9) and (10) gives for A}^: 

-exp{-i2l3x2). (16) 



ch{'yb) — i 




V y 77 \y ) "1 _ 


sh{^b) — exp {i [a{a + b) + j3a] } 


—ch{'~fb) + i 


(l:iy\ « _ ("i^ ^ 
\ y ) 1 \y) 7 


sh{'yb) + exp {i [a{a + b) — Pa]} 



In the same way for the KP-2 unit cell, from the corresponding nonhomogeneous system of 
equations (9) and (11) for the coefficients nj and Xjj we obtain 



cos{/3a) + i{l — y)^ — y^ sin{j3a) — exp[ia{a + b) + jb] 



cos{Pa) — i{l — y)^+y^ sin{Pa) + exp [ia{a + b) — 'yb] 



exp{-2^X2), (17) 



Ajj — 



ch{'yb) + i 






i) 


f] 


5/1(76) — exp{—i [a{a + b) + Pa]} 


—0/1(76) — i 








(L)i 

\y J 7 


sh{'jb) + exp{—i [a(a + 6) — Pa]} 



exp{-i2Px2). (18) 



Hence, for the both types of unit cells (7) and (8) it is seen from (15) and (17) that the wave 
functions (13) in the "region I" are different: ip^ 7^ ilr^\ and from (16) and (18) we see that the 
wave functions (14) in the "region 11" are also different: V'// 7^ V'//- Consequently, the choice 
of matching points inside the unit cell through relations (7)-(ll) determines the form of the 
periodic parts (13) and (14) of Bloch solutions for the model under consideration. 

Let's consider the case y = 1, i.e. when the electron effective mass (2) is constant: mi = 
1712 = fno- In this case P and 7 in (6) become 



Po = \l^-^E, 7o 



2m, 



-{V-E) 



(19) 



and equation (12) is reduced to the well-known Kronig-Pcnney equation [11]. If for the KP-1 
model (Fig. la) we set a;2 = a in relation (7), and then, making use of (16) and (19), we find for 
the coefficient X^j = {Xjj) the following value: 



^ o chi-job) ~ i^shi'job) - exp {i [a{a + b) + Poo]} 



-ch{'job) - i^shijob) + exp {i [a{a + b)- Poo]} 



exp{—i2Poa), 



(20) 



which coincides with the corresponding "expression (17)" in Bloss' paper [2]. For the KP-2 
model (Fig.lb) if X2 = in (8) we have the classical KP-model [11]. In this case from (17) and 
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(19) and from (18) and (19) we find the following expressions for the corresponding coefficients 
/.f^(/.|)^^andA?,^(Af,)^^: 

2\KP _ (^os{Poa) - ^sin{(3oa) - exp[ia{a + b) + -job] 

~ -cos{Poa) - ^sin{Poa) + exp[ia{a + b) - -job] ^ ' 

1.2 ,KP _ + '4^'''^^"^') - c:rp{-i\^^{a + h) + .7„fl]} 

—ch{'job) + i^sh{'job) + exp{—i[a{a + b) — f3oa\} 

These values (21) of (fij)^^ and (22) of (A|j)^^ have been implicitly used in [11]. However, 
choosing X2 = \b in (8), from (17) and (19) and from (18) and (19) we find that the corresponding 
coefficients )u| = (/Xj)*^ and A|j = (Af^)*^, given by 

2nG _ (^^iM - ^sin{Poa) - exp[ia{a + b)+ -fob] 

" -cosipoa) - fjiniPoa) + exp[ia{a + 6) - job] ""^^^"^"^^ ^^^^ 

^G ch{job) + i^sh{job) -exp{-i[a{a + b)+ (3oa]} 

Wi) = ^ exp{-tPob) (24) 

-c^job) + i^s^job) + exp {-i[a{a + b) - Poo]} 

are identical to Gubanov's results [12]. Comparing (21) with (23), and (22) with (24), we 
immediately obtain the following relations between them: 

if^^f = ifJ^jf^'expi-job), {Xhf = {\]i)''''exp{-iPob). 

Note that the exponential factors in (23) and (24) is due only to the different choice of the 
coordinate origin in the unit cell. 

The case of the Kronig-Penncy models with rectangular potential wells and partially constant 
effective mass (2) we treat in the same way, and the corresponding results are obtained taking 
into account the periodic potential V{x) (Fig.2): 

^ fo, Q}{x) -.nia + b) <x <nh+{n+l)a 
Vix) = < ^ (25) 

\-V, n^{x):nb+{n+l)a<x<{n+l){a + b), (n = 0, ±1, ±2, . . . ) 

where 1/ > 0. Relations (7) and (8) again specify two types of unit cells: the unit cell of the type 
"well-barrier" called KP-la type cell (Fig. 2a), and the unit cell of the type " barrier- well" called 
KP-2a type cell (Fig.2b). In the regions with constant potentials the quantities 9 and (/p, which 
now play the roles of 7 and (3 in (4) and (5), are related to the electron energy E by means of 
the equations 

In the case 0^ > 0, (^^ > the Schrodinger equation (3) for an electron with effective mass (2), 
moving in the potential field (25), for the both types of unit cells has nontrivial solutions only 
if the quantities 9 and satisfy the corresponding generalized Kronig-Penney relation 

COS [a(o + b)] = cos{9a)cos{(pb) + ^ " ^ ^ sin{9a)sin{<fb), (26) 

2y9(p 
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where 2/ = Now the periodic part of the Bloch wave function is of the form 

(jy{{x) = B ^^i^{exp[i{-a + e)x] + exp[i{-a - e)x]^ , (j = l,2), (27) 

in "regions I" (barriers): 0^(x) : n(a + b) < x < nb+ {n+ l)a, and of the form 

^^(x) = D |xijexp[z(-Q; + ip)x] + exp[i{-a - (p)x] | , (j = 1, 2) (28) 

in "regions IF (wells): 0^(x) : nb + (n + l)a < x < {n + l)(a + 5). For the KP-la unit cell, 
determined by condition (7), the corresponding nonhomogeneous system of equations (9) and 
(10) gives the following expressions for i/j and Xii' 



cos{(pb) — i 


(1- 




sin{iph) — exp {—i [a{a + b) + I 


M) 


—cos{ipb) + i 


a 




sin{ipb) + exp {—i [a{a + b) — 


ea]} 



exp{-i2ex2), (29) 



x]i 



cos{9a) — i 






sin{6a) — exp{i [a{a + b) + ipb]} 


—cos{9a) + i 


[\ y ) e \y) e 


sin{9a) + exp {i [a{a + 6) — ipb]} 



exp{—i2(px2) , (30) 



For the KP-2a unit cell, determined by (8), the coefficients z/| and Xii are obtained from the 
nonhomogeneous system of equations (9) and (11) are given by: 



cos{ipb) + i 


(1- 


if y ip 


sin{(fb) — 


exp{i [a{a + 6) + 9a]} 


—cos{(pb) — i 


a- 




sin{ipb) - 


- exp {i [a{a + 6) — 9a\} 



exp{-i29x2), (31) 



Xn 



cos{9a) + i 


M 




i) 


e 


sin{9a) — exp{—i [a{a + b) + ipb]} 


—cos{9a) — i 


[{¥. 


)f- 


(i: 


)« 


sin{9a) + exp{—i [a{a + b) — ipb]} 



exp{—i2ipx2)- (32) 



It is clearly seen that the corresponding wave functions (27) and (28) satisfy the inequalities 
4>] ^ cf)] and / 4)] I- 

If the electron energy is negative: E = — e, £ > and —V < E < 0, we have 9'^ = < 

and (f'^ = ^^{E + V) > 0, and denoting = k^, the generalized Kronig- Penney relation (26) 
becomes 



cos [a{a + b)] = ch{ka)cos{ipb) + 



(1 - yfa' - jy'k' + ip') 
2yipk 



sh{ka)sin{ipb). 



(33) 



In this case the quantity 7) (29) for the KP-la model becomes 



cos{pb) — 






sin{ipb) — exp [- 


-ia{a + 6) ± ka] 


—cos{ipb) + 


i(l 


ip ~^ y <fi 


sin{ipb) + exp 


—ia{a + b) ^ ka] 



exp{±2kx2), (34) 
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and from (30) we obtain also (xn) ^ = [Xii) = [xii] ^ 



ch{ka) — i 






1) 


k 


sh{ka) — exp {i [a{a + b) + ipb]} 


—ch{ka) + i 


[(¥ 


)i- 


il. 


H 


sh{ka) + exp{i [a{a + b) — ipb]} 



exp{—i2ipx2) ■ (35) 



For the KP-2a model the quantity (31) becomes 
(-1)" = 



cos{ipb) + 




sin{<pb) — 


exp [ia{a + 6) T ka\ 


—cos{ipb) — 


'iii-y)^Tyi 


sin{(pb) - 


- exp [ia{a + 6) ± ka] 



exp{±2kx2), 



(36) 



and from (32) we again have {xji)'^ = {xji) = {xji)^- 

ixhf = 



ch{ka) + i 


(¥)f + 




sh{ka) — exp {- 


i[a{a + b) + (pb]} 


—ch{ka) — i 


\ y ) k 


{y) ^ 


sh{ka) + exp {- 


-i [a{a + b) — fb]} 



exp{—i2ipx2) ■ 

(37) 



Hence, we have the relations 



When y = 1 our generalized Kronig-Penney relations (26) and (33) reduce to the corresponding 
Kronig- Penney relations given in [18]. 

In conclusion we can say that the differences between the wave functions (13) and (14) and 
between the wave functions (27) and (28), which finally depend on the choice of unit cell and 
the matching conditions, should be taken into account when we compute surface and interface 
electronic states. In other words, varying the position of the barrier-well interface within an unit 
cell, we will influence the conditions for the existence of surface and interface states. For the two 
potentials (1) and (25) our generalized Kronig-Penney relations (12) and (26), and the particular 
case (33), are different, as it is in the original Kronig-Penney relations, and the corresponding 
Bloch wave functions for these two potentials are also different. 
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Fig.l: One-dimensional Kronig-Penney model with rectangular barriers and lattice constant 
(a + 6): 6-tlie width of the barriers, T^-the height of the barriers: (a) unit cell of KP-1 type: 
"region II- region I" (well-barrier); (b) unit cell of KP-2 type: "region I- region 11" (barrier- well). 

Fig.2: One-dimensional Kronig-Penney model with rectangular wells and lattice constant 
(a -|- b): 6-thc width of the wells, IZ-the depth of the wells: (a) unit cell of KP-la type: "region 
II- region I" (well-barrier); (b) unit cell of KP-2a type: "region I- region 11" (barrier- well). 
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